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Abstract 

The quantum tunneling of the magnetization vector are studied theoretically 
in single-domain ferromagnetic nanoparticles placed in an external magnetic 
field at an arbitrarily directed angle in the ZX plane. We consider the mag- 
netocrystalline anisotropy with trigonal and hexagonal crystal symmetry, re- 
spectively. By applying the instanton technique in the spin-coherent-state 
path-integral representation, we calculate the tunnel splittings, the tunneling 
rates and the crossover temperatures in the low barrier limit for different an- 
gle ranges of the external magnetic field (Oh = tt/2, tt/2 <C Oh <C vr, and 
Oh = 7r). Our results show that the tunnel splittings, the tunneling rates and 
the crossover temperatures depend on the orientation of the external mag- 
netic field distinctly, which provides a possible experimental test for magnetic 
quantum tunneling in nanometer-scale single-domain ferromagnets. 
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I. INTRODUCTION 



Recently, there has been great experimental and theoretical effort to observe and interpret 
macroscopic quantum tunneling (MQT) and coherence (MQC) in nanometer-scale magnets 
at sufficiently low temperature.0 Theoretical investigations based on the spin-coherent-state 
path integral were performed for the single-domain ferromagnetic (FM) nanoparticles, which 
showed that MQT and MQC were possible in magnets containing as much as 10 5 — 10 6 
spins. Several experiments involving resonance measurements, magnetic relaxation, and 
hysteresis loop study for various systems showed either temperature-independent relaxation 
phenomena or a well-defined resonance depending exponentially on the number of total 
spins, which supported the idea of magnetic quantum tunneling.0 

More recently, the tunneling behaviors of the magnetization vector were studied exten- 
sively for the single-domain FM nanoparticles in the presence of an external magnetic field 
applied at an arbitrary angle. The MQT problem for FM particles with uniaxial crystal 
symmetry was first studied by Zaslavskii who calculated the tunneling exponent, the pre- 
exponential factors and their temperature dependences in the low barrier limit with the 
help of mapping the spin system onto a one-dimensional particle system.! For the same 
crystal symmetry, Miguel and Chudnovskyi calculated the tunneling rate by applying the 
imaginary-time path integral, and demonstrated that the angular and field dependences of 
the tunneling exponent obtained by Zaslavskii's method and by the path-integral method 
coincide precisely. They also discussed the tunneling rate at finite temperature and sug- 
gested experimental procedures.! Kim and Hwang performed a calculation based on the 
instanton technique for FM particles with biaxial and tetragonal crystal symmetryEI and 
Kim extended the tunneling rate for biaxial crystal symmetry to a finite temperature.il 
The quantum-classical transition of the escape rate for FM particles with uniaxial crystal 
symmetry in an arbitrarily directed field was investigated by Garanin, Hidalgo and Chud- 
novsky with the help of mapping onto a particle moving in a double-well potential.! The 
switching field measurement was carried out on single-domain FM nanoparticles of Barium 
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ferrite (BaFeCoTiO) containing about 10 5 — 10 6 spins.0 The measured angular dependance 
of the crossover temperature was found to be in excellent agreement with the theoretical 
predictionji which strongly suggests the MQT of magnetization in the BaFeCoTiO nanopar- 
ticles. Lii et al. studied the MQT and MQC of the Neel vector in single-domain antiferro- 
magnetic (AFM) nanoparticles with biaxial, tetragonal, and hexagonal crystal symmetry in 
an arbitrarily directed field.! 

In this paper, we extend the previous theoretical results obtained for the single-domain 
FM particles with biaxial and tetragonal symmetry to those for FM particles with a much 
more complex structure placed in an external magnetic field at an arbitrarily directed angle 
in the ZX plane, based on the instanton technique in the spin-coherent-state path-integral 
representation. We consider the magnetocrystalline anisotropies with trigonal and hexagonal 
crystal symmetry, respectively. Both the Wentzel-Kramers-Brillouin (WKB) exponents and 
the preexponential factors are evaluated analytically in the tunneling rates for MQT and the 
tunnel splittings for MQC in FM particles for different angle ranges of the external magnetic 
field {Oh = tt/2, vt/2 + O (e 3/2 ) < 9 H < n - 0(e 3/2 ), and 9 H = n), and the temperature 
which corresponds to the crossover from the thermal to the quantum regime is clearly shown 
for each case. Our results show that the distinct angular dependence, together with the 
dependence of the WKB tunneling rate and the crossover temperature on the strength of 
the external magnetic field, may provide an independent experimental test for the magnetic 
tunneling in single-domain FM nanoparticles. The calculations performed in this paper are 
semiclassical in nature, i.e., valid for large spins and in the continuum limit. We analyze the 
validity of the semiclassical approximation, and find that the semiclassical approximation is 
rather good for the typical values of parameters for single-domain FM nanoparticles. 

This paper is structured in the following way. In Sec. II, we briefly review the basic ideas 
of the MQT and MQC in single-domain FM particles. In Sees. Ill and IV, we study the 
quantum tunneling of the magnetization vector for FM particles with trigonal and hexagonal 
crystal symmetry in the presence of an external magnetic field applied in the ZX plane with 
a range of angles tt/2 < 6h < 7r. The conclusions are presented in Sec. V. In Appendix 
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A, we explain briefly the computation of the preexponential factors in the WKB tunneling 
rate, and then apply this approach to obtain the tunnel splittings for FM particles with 
trigonal crystal symmetry in a magnetic field applied perpendicular to the anisotropy axis 
(Oh = 7r/2) in detail. 



II. MQT AND MQC OF THE MAGNETIZATION VECTOR IN FM PARTICLES 

In this section we briefly review some basic ideas of MQT and MQC of the magnetization 
vector in single-domain FM nanoparticles, based on the instanton technique in the spin- 
coherent-state path integral 

The system of interest is a nanometer-scale single-domain ferromagnet at a temperature 
well below its anisotropy gap. For such a FM particle, the tunnel splitting for MQC or the 
tunneling rate for MQT is determined by the imaginary-time transition amplitude from an 
initial state \i) to a final state |/) as 

Ufi = (/I e~ HT \i) = J DQexp (S E ) , (1) 

where Se is the Euclidean action and DQ is the measurement of the path integral. In the 
spin-coherent-state path-integral representation, the Euclidean action can be expressed as 



(2) 

Tr = 



7 \dr J 7 \^ r , 

where V is the volume of the FM particle and 7 is the gyromagnetic ratio. M 
tvyS/V, where S is the total spin of FM particles. It is noted that the first two terms in Eq. 
(2) define the topological Berry or Wess-Zumino, Chern-Simons term which arises from the 
nonorthogonality of spin coherent states. The Wess-Zumino term has a simple topological 
interpretation. For a closed path, this term equals — iS times the area swept out on the unit 
sphere between the path and the north pole. The first term in Eq. (2) is a total imaginary- 
time derivative, which has no effect on the classical equations of motion, but it is crucial for 
the spin-parity effectslE However, for the closed instanton or bounce trajectory described 
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in this paper (as shown in the following), this time derivative gives a zero contribution to 
the path integral, and therefore can be omitted. 

In discussing macroscopic quantum phenomena, it is essential to distinguish between 
two types of processes: MQC (i.e., coherent tunneling) and MQT (i.e., incoherent tunnel- 
ing). In the case of MQC, the system in question performs coherent NH 3 -type oscillations 
between two degenerate wells separated by a classically impenetrable barrier. Tunneling 
between neighboring degenerate vacua can be described by the instanton configuration with 
nonzero topological charge and leads to a level splitting of the ground states.0 The tun- 
neling removes the degeneracy of the original ground states, and the true ground state is 
a superposition of the previous degenerate ground states. For the case of MQT, the sys- 
tem escapes from a metastable potential well into a continuum by quantum tunneling at 
sufficiently low temperatures, and the tunneling results in an imaginary part of the energy 
which is dominated by the so-called bounce configuration with zero topological charge. § As 
emphasized by Leggett, the two phenomena of MQC and MQT are physically very different, 
particularly from the viewpoint of experimental feasibility.^ MQC is a far more delicate 
phenomenon than MQT, as it is much more easily destroyed by an environment,^ and by 
very small c-number symmetry breaking fields that spoil the degeneracy. 

In the semiclassical limit, the dominant contribution to the transition amplitude comes 
from the finite action solution (instanton) of the classical equation of motion. The motion 
of the magnetization vector ~R% is determined by the Landau-Lifshitz equation, 



which can also be expressed as the following equations in the spherical coordinate system, 



where 9 and <p denote the classical path. Note that the Euclidean action Eq. (2) describes 
the (1 © l)-dimensional dynamics in the Hamiltonian formulation with canonical variables 
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(3) 




(4b) 



(4a) 



and Pcf, = S(l— cos#). The instanton's contribution to the tunneling rate T for MQT or 
the tunnel splitting A for MQC (not including the topological Wess-Zumino or Berry phase) 



where ui p is the frequency of small oscillations near the bottom of the inverted potential, and 
S c i is the classical action. The preexponential factor A originates from the quantum fluctu- 
ations about the classical path, which can be evaluated by expanding the Euclidean action 
to second order in the small fiuctuations.0'0 In Ref. 12, Garg and Kim studied the general 
formalism for calculating both the exponent and the preexponential factors in the WKB 
tunneling rates for MQT and MQC in single-domain FM nanoparticles. In Appendix A, we 
explain briefly the basic idea of this calculation, and then apply this approach to calculate 
the instanton's contribution to the tunnel splittings for MQC of the magnetization vector 
in FM particles with trigonal crystal symmetry in an external magnetic field perpendicular 
to the anisotropy axis (considered in Sec. Ill) in detail. 

III. MQC AND MQT FOR TRIGONAL CRYSTAL SYMMETRY 

In this section, we study the tunneling behaviors of the magnetization vector in single- 
domain FM nanoparticle with trigonal crystal symmetry. The external magnetic field is 
applied in the ZX plane, at an angle in the range of tt/2 < 9h < 7T. Now the total energy 
E (9, 0) can be written as 

E {6, 0) = Ki sin 2 9 - K 2 sin 3 9 cos (30) - M H X sin 9 cos - M H Z cos 9 + E , (6) 

where K\ and K 2 are the magnetic anisotropy constants satisfying K\ ^> K 2 > 0, and E 
is a constant which makes E (9, 0) zero at the initial orientation. As the magnetic field is 
applied in the ZX plane, H x = Hsm9n and H z = H cos 9h, where H is the magnitude of 
the field and 9h is the angle between the magnetic field and the z axis. 



is given 



byffl 




(5) 
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In the absence of the external magnetic field, the system reduces to one with threefold 
rotational symmetry around z axis and reflection symmetry in the XY plane. The unit vec- 
tors z and — z define the two classical ground state configurations. The transition amplitude 
between degenerate ground states can be suppressed to zero resulting from the destructive 
Wess-Zumino phase if the system has time-reversal invariance at zero magnetic field.0 How- 
ever, for the closed instanton or bounce trajectory described in this paper (as shown in the 
following) the phase term in Eq. (2), proportional to dip/dr (not (dtfi/dr) cosO term) gives 
a zero contribution to the integral Eq. (2) and, therefore, can be omitted. 

By introducing the dimensionless parameters as 

K 2 = K 2 /2K X ,H X = H x /H , H z = H z /H , (7) 

Eq. (6) can be rewritten as 

E (9, 0) = - sin 2 9 - K 2 sin 3 9 cos (30) - H x sin 9 cos - H z cos 9 + E Q , (8) 

where E (9, 0) = 2K\E (9, 0), and H = 2K\/M G . At finite magnetic field, the plane given 
by = is the easy plane, on which E (9, 0) reduces to 

E {6, = 0) = - sin 2 6 - K 2 sin 3 9-Hcos(9- 6 H ) + E . (9) 

We denote #o to be the initial angle and 9 C the critical angle at which the energy barrier van- 
ishes when the external magnetic field is close to the critical value H c (9h) (to be calculated 
in the following). Then, the initial angle 9q satisfies [dE (9, = 0) /d6\ e _ e = 0, the critical 
angle 9 C and the dimensionless critical field H c satisfy both \dE {9, = 0) /d9\ e _ Q = 

and [d 2 E(9,(p = 0) /^ 2 ] e=ec _ = _ = 0, which leads to 

- sin 

(20 o ) - 3K 2 sin 2 9 cos 9 + H sin (9q — 9h) = 0, (10a) 

- sin (29 c ) - 3K 2 sin 2 9 C cos 9 C + H c sin (6 C - 9 H ) = 0, (10b) 
cos (26 c ) - 3K 2 (2 sin 9 C cos 2 9 C - sin 3 9 C ) + H c cos (9 C - 9 H ) = 0. (10c) 

After some algebra, H c (9h) and 9 C are found to be 
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1 

(sin^) 2/3 + |cos^| 2/3 
+ 6 ^- " —2 



3/2 



1 + 3fT 2 - 



(l + |cot^| 2/3 ) 



1/2 



(l + |cot^| 2/3 ) 



sin 9 r 



1 



1 + Icotfln 



2/3 



1-2K 



^ I cot e H \ 2/3 



^ | cot 6 H \ 2/3 



3/2 



(11a) 



(lib) 



(l + |cot0 ff | 2/3 ) (l + |cot0 H | 2/3 )' 

Now we consider the limiting case that the external magnetic field is slightly lower than 
the critical field, i.e., e = 1 — H/H C <C 1. At this practically interesting situation, the barrier 
height is low and the width is narrow, and therefore the tunneling rate in MQT or the tunnel 
splitting in MQC is large. Introducing r\ = 9 c — 9 (\r)\ <C 1 in the limit of e 1), expanding 
[dE (9, = 0) /d9] e=6o = about 9 C , and using the relations [dE (9, = 0) /dQ\ e=e jj = -jj = 
and [d 2 E {9, = 0) /g?0 2 ] q _ q = 0, we obtain the approximation equation for 77 in the 

order of e 3 / 2 , 



-eH c sin (9 C - 9 H ) - rj 2 ( - sin 29 c + 3K 2 cos 39, 



+V 



eH c cos (9 C -9 H )+rf [- cos 29 c - 3K 2 sin 39, 



0. 



(12) 



Then E (9, 0) reduces to the following equation in the limit of small e, 



E (5, 0) = 2K 2 sin 2 (30/2) sin 3 (0 O + S) + H X sin (0 O + 5) (1 - cos 0) + ^ (5) , (13) 
where 5 = 9 — 9 (\5\ <C 1 in the limit of e <C 1), and E 1 (5) is a function of only 5 given by 



Ei (6) = 



if c sin (9 C — 9 H ) — K 2 ( cos 3 9 C — - sin 2 9 C cos 0, 



1 



— \H C cos (0 C - 9 H ) - 3fT 2 (sin 3 C - 4 sin C cos 2 C ) ] 

3 

--X 2 (sin 3 c -4sin0 c cos 2 c ) 5 2 e. 



(3S 2 V - 5 s ) 



(14) 



In the following, we will investigate the tunneling behaviors of the magnetization vector 
in FM particles with trigonal crystal symmetry at different angle ranges of the external 
magnetic field as 9h = vr/2 and n/2 < 9h < 7r, respectively. 



A. 6 H = tt/2 

For 6 H = tt/2, we have 6 C = tt/2 from Eq. (lib) and rj = ^2e (l + 9K 2 /2) from Eq. 
(12). Then E l (5) of Eq. (14) reduces to 

Ei (8) = (5 - 2 V ) 2 . (15) 

The plot of the effective potential E x (5) as a function of 5 (= 9 — 6 ) ioi9 H = tt/2 is shown in 
Fig. 1. Now the problem is one of MQC, where the magnetization vector resonates coherently 
between the energetically degenerate easy directions at 5 = and 5 = 2y / 2e (l + 9K 2 /2) 
separated by a classically impenetrable barrier at 5 = \[2e (l + 9K 2 /2). Substituting Eq. 
(15) into the classical equations of motion, we obtain the classical solution called instanton 
as 

5 = V2~e (l + p( 2 ^J [1 + tanh (cJ c r)} , (16) 

where uJ c = a/c/2 (l + 21K 2 /2 — e/2), r = uj t, and uj = 2KiV/hS. We can calculate 
the classical action by integrating the Euclidean action Eq. (2) with the above classical 
trajectory, and the result is found to be 

S d = (l + ^K 2 + ^ . (17) 

Now we consider the transition exponent which is usually addressed by experiments. 
Transitions between two states in a bistable system or escaping from a metastable state can 
occur either due to the quantum tunneling or via the classical thermal activation. In the 
limit of temperature T — > 0, the transitions are purely quantum-mechanical and the rate 
goes as T ~ exp (—Sd), with Sd being the classical action or the WKB exponent which is 
independent of temperature. As the temperature increases from zero, thermal effects enter 
in the quantum tunneling process. If the temperature is sufficiently high, the decay from a 
metastable state is determined by processes of thermal activation, and the transition rate 
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follows the Arrhenius law, T ~ exp (—U/ksT), with ks being the Boltzmann constant and 
U being the height of energy barrier between the two states. Because of the exponential 
dependence of the thermal rate on T, the temperature T c characterizing the crossover from 
quantum to thermal regime can be estimated as ksT c = U/S c i. For a quasiparticle with the 
effective mass M moving in one-dimensional potential U (x), a more accurate definition of 
the crossover temperature in the absence of any dissipation was presented by GoldanskiiJliE 
ksT'c = HuJb/^TT, where Ub = \/—U" (xb) /M is the frequency of small oscillations near the 
bottom of the inverted potential, —U (x), and Xb corresponds to the bottom of inverted po- 
tential. Below T' c , thermally assisted quantum tunneling occurs from the excited levels, that 
further reduces to the quantum tunneling from the ground-state level as the temperature 
decreases to zero. Above T' c , quantum tunneling effects are small and the transitions occur 
due to the thermal activation to the top of the barrier. For the MQT problem, i.e., the 
problem of decay from the metastable state, both T c and T' c can be used as the definition of 
the crossover temperature corresponding to the crossover from classical to quantum behavior 
since the quantum escaping from a metastable state is one process of incoherent tunneling. 
However, for the MQC problem, i.e., the problem of resonance between degenerate states, 
the situation is different. As the temperature growing from zero, three kinds of transitions 
should be taken into account: quantum coherence between the degenerate ground-state 
levels (coherent tunneling), quantum tunneling from the excited levels (thermally assisted 
tunneling or incoherent tunneling), and classical over-barrier transition (incoherent transi- 
tion). Two kinds of crossover temperatures can be defined to distinguish the three regimes. 
The Goldanskii definition T' c for MQC problem corresponds to the crossover from quantum 
coherence between the degenerate ground-state levels (coherent tunneling) to quantum tun- 
neling from the excited levels (thermally assisted tunneling or incoherent tunneling), while T c 
corresponds to the crossover from quantum coherence between the degenerate ground-state 
levels (coherent tunneling) to classical over-barrier transition (incoherent transition). Exper- 
iments involving magnetic relaxation and resonance measurements for various systems have 
shown either temperature-independent relaxation phenomena (in MQT) or a well-defined 
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resonance (in MQC) below some crossover temperature, which strongly support the exis- 
tence of quantum tunneling processes.0 And more recently, the crossover from quantum to 
classical behavior and associated phase transition have been investigated extensively in MQT 
and MQC in single-domain FM particles.0Hll It is noted that the sharpness of the crossover 
between thermal and quantum regimes also depends on the strength of the dissipation with 
environment. In the case of the low dissipation which is common for the magnetic systems, 
its effect on the crossover is smallliB 

For the single-domain FM nanoparticle in a magnetic field applied at Oh = tt/2, 
the magnetization vector resonates coherently between the energetically degenerate easy 
directions at 5 = and 5 = 2^/2e (l + 9K 2 /2) separated by a classically impenetra- 
ble barrier at 5 = \/2e (l + 9K 2 /2), and the height of energy barrier is found to be: 
U = KiVe 2 (l + 1%K 2 \ Then, equating S c i to U/ksT, we obtain that the crossover from 
quantum coherence between the degenerate ground-state levels (coherent tunneling) to clas- 
sical over-barrier transition (incoherent transition) occurs at 



, m 3 1/2 K l V ( 21— 1 \ 

k B T c = -^t 1 ' 2 ^— 1 H K 2 - -e . 

2 5/2 g y 2 2 2 J 



(18) 



For this MQC problem, the Goldanskii definition T' c corresponding the crossover from quan- 
tum coherence between the degenerate ground-state levels (coherent tunneling) to quantum 
tunneling from excited levels (thermally assisted tunneling or incoherent tunneling) be- 



comes ksT'c = HuJb/^TT, where Ub = u^b^o, with ZJb = y — E (5 m ) /M is the frequency of 
small oscillations of the magnetization vector near the bottom of the inverted potential, 
M~ l = (l + 12K 2 — e), and 5 m is the position of the energy barrier. For the present case, 
5 m = V2e (l + 9K 2 /2) and U b = V^( 1 + 21^ 2 /2 - e/2) = V2ZJ C - Then it is easy to obtain 
that 

(*+!*-?)■ 

The comparison of Eqs. (18) and (19) shows that T c w 1.67T C ', which is consistent with the 
physical interpretation for quantum-classical transition in the MQC problem. 
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It is noted that the quantum tunneling of the magnetization vector in single-domain 
FM nanoparticles are studied with the help of the instanton technique in the spin-coherent- 
state path- integral representation, which is semiclassical in nature, i.e., valid for large spins 
and in the continuum limit. Therefore, one should analyze the validity of the semiclassical 
approximation. It is well known that for this approach to be valid, the tunneling rate must be 
small, which indicates that the WKB exponent or the classical action Sd 3> 1. Moreover, the 
energy huf, of zero-point oscillations around the minimum of the inverted potential — ~E\ (5) 
should be sufficiently small compared to the height of the barrier, U = 2K{VEi (8 m ) . 
For the single-domain FM nanoparticle with trigonal crystal symmetry in a magnetic field 
applied at 9h = tt/2, it is easy to show that the WKB exponent is approximately given by 



which agrees up to the numerical factor with the result of the classical action in Eq. (17) 
obtained by applying the explicit instanton solution. For the typical values of parameters 
for single-domain FM nanoparticles, K\ ^ 10 8 erg/cm 3 , K 2 ^ 10 5 erg/cm 3 , and the total 
spin S = 10 6 , we obtain that B - U/ftw b « 15.8 from Eq. (20) and S d « 59.6 for e = 0.001 
from Eq. (17). In this case the semiclassical approximation should be already rather good. 

By applying the instanton technique for FM particles in the spin-coherent-state path- 
integral representation, we obtain the instanton's contribution to the tunnel splitting as 
(for detailed calculation see Appendix A), 



where the WKB exponent or the classical action S c i has been presented in Eq. (17). 

Now we apply the effective Hamiltonian approach to evaluate the ground-state tunnel 
splitting.El For the present case, the effective Hamiltonian can be written as 




(20) 









—HA, 



o 



(22) 



-ft An 
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A simple diagonalization of H e ff shows that the eigenvalues of this system are ±hA . There- 
fore, the splitting of ground state due to resonant coherently quantum tunneling of the mag- 
netization vector between energetically degenerate states is HA = 2hA , where hA is shown 
in Eq. (21) with Eq. (17) for single-domain FM particles with trigonal crystal symmetry in 
a magnetic field applied perpendicular to the anisotropy axis (9h = tt/2). 

B. vr/2 <6 h <tt 



For 7r/2 < 9 H < 7r, the critical angle 6 C is in the range of < 9 C < n/2, and n a/2c/3. 
Then E l {5) of Eq. (14) reduces to 

1 I cot 6 H \ 1/3 



E l (S) 



21 + 1 cot e H \ 2/3 



2 2 



1 

|cotfl~' 2/3 



1/2 



6e5 z - 5' 



(23) 



1 -r r< >i (/// 

The dependence of the effective potential E\ (5) on 5 (= 6 — 9 ) for 6h = Sir/ 4 is plotted in 
Fig. 2. Here, K 2 = 0.001. Now the problem becomes one of MQT, where the magnetization 
vector escapes from the metastable state at 5 = 0, = through the barrier by quantum 
tunneling. Substituting Eq. (23) into the classical equations of motion, the classical solution 
called bounce is found to be 



= i (6e) 3/4 | cot 9 H \ 1/6 (l + | cot 9 H \ 2/3 J 



x 1/2 



l + |-|^ 2 (l + |cot^| 2 / 3 ) 



1/2 



, K 2 2|cot^| 2/3 -9 - 
I — ; — : — — + X 2 - 



Icot^wl 273 -3 



(l + |cot^| 2/3 ) 
5 = V6e/ cosh 2 (cJ c r) , 



1/2 



1 + Icotfljy 



2/ 3 \ 



3/2 



sinh (u c t) 
cosh 3 (uJ r r) ' 



(24) 



which corresponds to the variation of S from S = at r = — ooto the turning point S 
at t = 0, and then back to 5 = at r = +oo, where 

ZJ r = 3 1 / 4 x 2" 3 /V/ 4 - 



6e 



I cot H 1 1/6 



1 + | cot H \ 2/3 

1 2/3 



l- 6 - + 9 -K 2 (\ f kotfl--' 2 " 



\ 1/2 



ir 2 2|cot^| z/d -21 _ |cot^| 2/3 + 3 



+ ^2 



(l + |cot^| 2/3 ) 1/2 (l + |cot^| 2/3 ) 



3/2 
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The associated classical action is then given by 



g d = 31/4x217/ V / 4 icot^r /6 



K 2 |cotfl H | 2/3 + 9/2 
171 



l + |-^ 2 (l + |cot^| 2/3 ) 



1/2 



I COt0 



1 2/3 



2 (l + |cot^| 2/3 ) (l + |cot^| 2/3 ) 

For this case, the barrier height is 



3/2 



(25) 



U = 2K l VE l (5 m ) 

_2 7 / 2 I cot 9 H \ 1/3 
= 1 + , cot ^ ,2/3 



(l + |cot ^| 2/3 ) 



1/2 



6 3/2 (^V) , 



at 5 m = 2(6e) 1 ^ 2 /3, and the frequency of small oscillations of the magnetization vector 



around the bottom of the metastable well is 

ZJ b = 3 1 / 4 X 2 l/4 g l/4 l COt ^! 1/6 

1 + |cot0 H | 2/3 



e 9— / l2/ o\ 1/2 

1-2 + 2^ (l + |cot^| 2/3 ) 



K 2 2|cot^| 2/3 -21 — 
+ — ' + #2 

(l + |cot^| 2/3V 

2uL. 



I cot 6 H \ 2/3 + 3 



1 + Icot 



2/3 



3/2 



Then the WKB exponent or the classical action B is approximately given by 
U 



B 



fai! b 
2 9/4 



5e 5/4 |cot 6 H 



1/6 



1 + \ -^2(1 + M e H \ 2 ' 3 ) 



37/4 

^2 I cot 0h 1 2/3 + 9/2 



1/2 



Icot 0h| 2/3 + 3 



(26) 



(1 + |cot e H \ 2 '*) 1/2 "(i + |cot^| 2 / 3 ) 3/2 

which is consistent with Eq. (25) up to the numerical factor. After a simple calculation, we 
obtain the crossover temperature as 



k B T c 



1/A K 1 V \cot6 H 



1/6 



23/4x3^/4 S l + |cot^| 2/3 
u_z — + K 2 



1 - 1+1*2(1+ icot ^r) 



1/2 



KT 2 |cot^| 2/3 -21/2 |cot^| 2/3 + 3 



(l + |cot^| 2/3 ) 1/2 (l + |cot^| 2/3 ) 



3/2 



(27) 
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corresponding to the transition from quantum to thermal regime. For a nanometer-scale 
single-domain FM particle, the typical values of parameters for the magnetic anisotropy 
coefficients are K\ = 10 s erg/cm 3 , and K 2 = 10 5 erg/cm 3 . The radius of the FM particle 
is about 12 nm and the sublattice spin is 10 6 . If e = 0.001, we obtain that T c (135°) ^ 
203mK corresponding to the crossover from quantum to classical regime, which compares 
well with the experimental result of 0.31K on single-domain FM nanoparticles of Barium 
ferrite (BaFeCoTiO).i Note that, even for e as small as 10~ 3 , the angle corresponding to an 
appreciable change of the orientation of the magnetization vector by quantum tunneling is 
S 2 = V6e rad> 4°. 

The classical action S c i can be obtained by solving numerically the equations of motion 
(4a) and (4b). In Fig. 3 we present the Oh dependence of S c i with e = 0.001 and K 2 = 0.001 
for 7r/2 < Oh < vr by numerical and analytical calculations, respectively. As is noted in the 
figure, the analytical result obtained from Eq. (25) is almost valid in the whole range of 
angles n/2 < H < 7r. 

By applying the formulas in Ref. 12, and using Eq. (25) for the WKB exponent or the 
classical action, we obtain the tunneling rate T of the magnetization vector in single-domain 
FM nanoparticles with trigonal crystal symmetry in a magnetic field applied in the range of 
7r/2 < Oh < 7T as 



i-i + |7? 2 (i + |cote H |^ 1/2 

e" Sci . (28) 



* 1/2 h 1 l + |cot^| 2/3 

K 2 |cot^| 2/3 -51/2 K 2 |cot^| 2/3 + 3 



(1+ icot^ 3 ) 172 2 (i+ icot^r) 3/2 



IV. MQC AND MQT FOR HEXAGONAL CRYSTAL SYMMETRY 

In this section, we study the quantum tunneling of the magnetization vector in single- 
domain FM particles with hexagonal crystal symmetry whose magnetocrystalline anisotropy 
energy E a (0, 0) at zero magnetic field can be written as 
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E a (0, 0) = K x sin 2 + K 2 sin 4 + K 3 sin 6 - K' 3 sin 6 cos (60) , (29) 

where Ki, K 2 , K 3 , and K' 3 are the magnetic anisotropic coefficients. The easy axes are ±z 
for K\ > 0. When we apply an external magnetic field at an arbitrarily directed angle in 
the ZX plane, the total energy of this system is given by 

E (0, 0) = E a (0, 0) - M H X sin cos - M H Z cos 9 + E , (30) 

By choosing K 3 > 0, we take = to be the easy plane, at which the potential energy can 
be written in terms of the dimensionless parameters as 

E(9, = 0) = - sin 2 9 + K 2 sin 4 9 + (k 3 - lf 3 ) sin 6 - tfcos (0 - 9 H ) + E , (31) 

where K 3 = K 3 j2K x and K 3 = K' 3 /2K 1 . 

Then the initial angle 0o is determined by [dE(9,0) /d9\ e _ e ^ = 0, and the critical an- 
gle 9 C and the dimensionless critical field H c by both [dE (9, 0) /d9] 0=0 jj = jj = and 
[d 2 E (9, 0) /d9 2 ] 0=dclI=T j c = 0, which leads to 

^ sin (20 o ) + # sin (9 - 9 H ) + 4K 2 sin 4 O + 6 (Z 3 - lf 3 ) sin 5 O cos 9 = 0, (32a) 

- sin (20 c ) + # c sin (0 C - H ) + 4Z 2 sin 4 C + 6 3 - lf 3 ) sin 5 9 C cos C = 0, (32b) 
cos (29 c ) + H c cos (0 C - H ) + AK 2 (3 sin 2 C cos 2 9 C - sin 4 C ) 

+6 (K 3 - K'^j (5 sin 4 9 C cos 2 C - sin 6 9 C ) = 0, (32c) 



Under the assumption that |-K"2 , 
if c as 

1 



if a - K, 



< 1, we obtain the dimensionless critical field 



H r = 



(sin0^) 2/3 + |cos0^| 2/3 



3/2 



1 + 



6 K 3 - K, 



1 + |cot0H 



2/3 



(l + |cot^| 2/3 )' 



(33) 



In the limit of small e = 1 — H/H c , Eq. (32a) becomes 

-e~H c sin (0 C - H ) + rf [(3/2) # c sin (9 C - 9 H ) + 3K 2 sin (40 c ) 
+12 3 - Z 3 ) sin 3 C cos C (5 - 8 sin 2 C ) + rj {eH c cos (0 C - 9 H ) 
-T] 2 [(1/2) # c cos (0 C - H ) + 4Z 2 cos (40 c ) 
+12 (F 3 - K'^j sin 2 C (5 - 20 sin 2 9 C + 16 sin 4 C )] } = 0, 
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(34) 



where r] = C — O which is small for e 1. By introducing a small variable 5 = — 
(\5\ <C 1 in the limit of e <C 1), the total energy becomes 

E (5, 0) = K 3 [1 - cos (60)] sin 6 (0 O + S) + H x (1 - cos 0) sin (0 O + 5) + ^ (5) , (35) 

where £i (<5) is a function of only S given by 



£i (5) 



l -H c sin (0 C - H ) + ^2 sin (40 c ) + a(k 3 - K 3 ^j (5 sin 3 C cos 3 C - 3 sin 5 C cos C ) 

x (<5 3 - 35 2 r]) + cos (0 C - H ) + Z 2 cos (40 c ) + 3 (Z 3 - lf 3 ) sin 2 C (sin 4 C 

-10 sin 2 C cos 2 C + 5 cos 4 C )] (<5 4 - 4<5 3 r/ + 65V - 45 2 e) + e<5 2 [4X 2 cos (40 c ) 



+ 12 ^ 3 - K'^j sin 2 C (sin 4 C - 10 sin 2 C cos 2 C + 5 cos 4 C ) 



(36) 



In the following we investigate the MQC and MQT of the magnetization vector in FM 
particles with hexagonal crystal symmetry for different angle ranges of the external magnetic 
field: H = tt/2, tt/2 + O (e 3 / 2 ) < 9 H < n - O (e 3 / 2 ), and 6 H = tt, respectively 



A. 6 H = vr/2 

For 9h = tt/2, i.e., the external magnetic field is applied perpendicular to the anisotropy 
axis, we obtain that C = n/2 and i] = \f2t |\ — 4K 2 — 12 ^3 — K 3 ^j . Now E\ (5) becomes 

E x (5) = ^ [l + 12Z 2 + 30 (X 3 - ^3)] ^ 2 {5 - Sv 7 ^ [l - AK 2 - 12 (Z 3 - Z 3 )] }' . (37) 

Substituting Eq. (37) into the classical equations of motion, we obtain the following instan- 
ton solution 



= ie 
5 = 



1 + e - - AK 2 - 18K 3 -6(k 3 - K' 3 ) 1 1 

2 V /J cosh (uj c r) 

1-AK 2 - 12 (K 3 - F' 3 )] [1 + tanh (cJ c r)] , 



(38) 



which corresponds to the variation of 5 from 5 = at r = —00 to 5 



2V2e 



1 - AK 2 - 12 [K 3 - K 3 







at t = +00, where 
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1--+AK 2 + 18K 3 + 6[K 3 -K 3 



We can calculate the classical action by integrating the Euclidean action of Eq. (2) with the 
above instanton solution, and the result is found to be 

2 5/2 



Sr. 



_ Se S/2 



1 + - - 8K 2 - 18K 3 -24 if 



K, 



(39) 



From Eq. (37) we obtain that the height of barrier is U = 2K\VEi (5 n 



1 - AK 2 - 18 \K* - K 



at <L 



1 — 4Ko — 12 I K-i — K 



and the os- 



cillation frequency around the minimum of the inverted potential — E\ (5) is 



1- 1+4^ + 184 + 6 ( Is :, 



K' 



v2uj c . 



Then the WKB exponent is approximately given by 

U 



B 



1 



1 + - - 8K 2 - 18K 3 -24[K 3 -K 3 



(40) 



which agrees up to the numerical factor with Eq. (39) obtained by applying the explicit 
instanton solution. The temperature corresponding to the crossover from the quantum 
coherence between the degenerate ground-state levels (coherent tunneling) to the classical 
over-barrier transition (incoherent transition) is found to be 



1-- + 4K 2 + 18K 3 + 6[K 3 -K 3 



(41) 



2 5 / 2 S 

and the temperature corresponding to the crossover from quantum coherence between the 
degenerate ground-state levels (coherent tunneling) to quantum tunneling from excited levels 
(thermally assisted tunneling or incoherent tunneling) is found to be 



7T D 



1 - - + 4K 2 + 18K 3 + 6 ( K 3 - K 3 



(42) 



By applying the instanton technique for single-domain FM particles in the spin-coherent- 
state path-integral represent at ion, @ we obtain the instanton's contribution to the tunnel 
splitting, ftA as 
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2 13/4 



7T 



1/2 



l-- + 9K 3 -6[K 3 -K 3 



(43) 



where the WKB exponent or the classical action S c i is clearly shown in Eq. (39). Then the 
splitting of ground state due to resonant coherently quantum tunneling of the magnetization 
vector between energetically degenerate states is found to be HA = 2hA for FM particles 
with hexagonal crystal symmetry in a magnetic field applied perpendicular to the anisotropy 
axis (&h = 7r/2) with the help of the effective Hamiltonian approach. 



B. vr/2 + O (e 3 / 2 ) < 9 H < ir - O (e 3 / 2 ) 



For this case, r\ ps y2e/3 and the critical angle 9 C is found to be 



sin 6 C = 



(l + |cot^| 2/3 ) 
Now Ei (5) becomes 



1/2 



8 \coiO H \ 2/S " 2/3+ s(K 3 -K' 3 ) ' COt ^' 2/3 , 



Ei (5) = 



i | cot e H | 1/3 

2l + |cot^| 2/3 



2/3 



. 4_7-4|cot0 H 
1 - : -^Ta 



3 l + |cot0#| 

1 2/3 



c (~F7 T7>\ 11 - 16 Ic0t6 l ^| 2,/3 / r— r2 r o\ 

+2( < K *-^ ( 1 + l eo t ^f j 
Then the classical equations of motion have the following bounce solution 



= i (6e) 3/4 |cot^| 1/6 (l + |cot^| 2/3 ) 



1/2 



2/3 



1 e 4— 5- |cot0#, 
2 3 l + |cot^| 2/3 



-18^- 



' 1 + |cot 
5 = VQe/ cosh 2 (uj c t) 



i -2(k 3 -k' 3 ) '-^'-r 



sinh (u> c t) 
cosh 3 (cJ c r) ' 



(44) 



(45) 



corresponding to the variation of S from 5 = at r = — oo to the turning point <5 = y6e at 
r = 0, and then back to 5 = at r = +oo, where 
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'3\ 1/4 1/4 \cot6 H \ 1/6 [ e 4 5-3|cot^| 2/3 
6 l + |cot^| 2/3 [ 2 + 3 2 l + |cot^| 2/3 



+ I8X3- 



1 1 o/V ^ 7-10|cot^| 2 / 3 

i i + icot^r +2 ^ 3 -^ ( 1+ | COt ^/3) 2 

The classical action associated with this bounce solution is found to be 

_ 2^x3V 4 5/4 ,i/6 [ e 4 2-lcot^i 2 / 3 

5c/ " 5 * |COt ^' 1 + 2 + 3^ 2 1 + | cot ^/3 



-1SK' 3 - 



1 2/3 



1 / ,\ 2-3|cot0 £ 

l + |eot^| 2/3 V 7 (l + |cot^| 2/3 ) 

For this case, the barrier height C/ (= 2ii'il/£'i [5 m = 2\/6e/3)) is given by 

2 7/2 



u = 



l + |cot^| 2/3 



V ^ (l + |cot9„| 2 ' 3 ) 

and the frequency of small oscillations of the magnetization 
the inverted potential —E x is 



4 7-4|cot^| 2/3 
3 l + |cot^| 2/3 

2/3 



(46) 



vector around the minimum of 



- (5) is 

zj b = 3 1 / 4 x 2 v4 e i/4 \™te H \ 1/6 

1 + Icotflwl 



l + |cot^| 2/3 [ 2 3 



+ I8X3 

= 2ZJ C . 



, e 4 5-3|cot^| 2/3 
l + |cot^| 2/3 

2/3 



1 + |cot.e„|=/ 3 ^ ^(i + icotfei 2 ") 



Then the WKB exponent is approximately given by 
U 



B 



2 9/4 



Se 5 / 4 |cot0tf| 1/6 



e 4 _2-|cot0 H 
J- + - + -^2 



2/3 



-18X, 



2 ' 3^i + | co t^| 2/3 

2-3|cot ^| 2/3 



4 (X 3 - 



i + icot^i 2/3 1 ^ 3 li3y ri + icot^i 2 / 3 ) : 



(47) 
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which agrees with Eq. (46) up to the numerical factor. Equating the classical action S c i 
to U/kBT c , where U is the barrier height, we obtain that the crossover from quantum to 
classical behavior occurs at 

5 1/4K1V |cot0#| 1/6 



k B T c 



2 3 /4 x 3 7 /4 



+ 18K* 



S l + |cot0#| 2/3 



e 4—5-3 \cot6 H \ 2/3 

1 ~~ o + 0-^2 ; -777- 

2 3 1+ cot ^ 2/3 



1 + |cot#tf| 2/3 



|cot#~' 2/3 



(48) 



1 "t CO l Off 1 

Based on the instanton technique,0 we obtain the tunneling rate corresponding to the 
escaping of the magnetization vector from the metastable state for single-domain FM 
nanoparticles with hexagonal crystal symmetry in a magnetic field applied in the range 
of it/2 + O (e 3 / 2 ) < 6 H < 7r — O (e 3 ^ 2 ) as the following equation, 



2 3 V8 x 37/8 v KiS . 1/2eV8 \cot9 H 



1/4 



7T 



1/2 



h 



1 + Icot0ff 



2/3 



1 - 7 + 9K' 3 m 

4 3 l+cot0 H | 2/3 



2— 12 - 7|cote H 

A 

3 1 



2/3 



cot e r { 2/3 



7 H\ 



(l + \cot0 H \ 



2/3 



-S c l 



(49) 



where the WKB exponent or the classical action S c i has been clearly shown in Eq. (46). 



C. 9 H = 7T 

Finally, we study the MQT of the magnetization vector corresponding to the escaping 
from the metastable state in single-domain FM nanoparticles with hexagonal crystal sym- 
metry in a magnetic field applied at 6 H = it, i.e., antiparallel to the anisotropy axis. Now 
the total energy become 



E(5A) = K' 3 [l-cos(6(f>)]5 6 + ±5 2 
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1 - 32AT 2 + 48 (K 



K' 



6- 



(50) 
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The classical equations of motion have the bounce solution 



— U7T 

<f> = -IU C T + — , 



5 = 



\ 



4e 



lKo-6 



32K 3 (1 - cosh (6cJ c r)) + § (1 - A8K 2 + 96 [K 3 - K 



(51) 



where n — 0, 1, 2, 3, 4, 5, and u> c = e. The corresponding classical action is found to be 



Sri = -Se—— In \ 

3 Ax I A 2 



with 



(52) 



Ai = 1 - 8^ - 32fT 3 e - -e 

3 



1 - 48tf 2 + 96 (tf 3 - , 



(53) 



and 



A 9 = 32K,e. 



(54) 



According to the formulas in Ref. 12, we obtain the tunneling rate of the magnetization vec- 
tor escaping from the metastable state for single-domain FM nanoparticles with hexagonal 
crystal symmetry in a magnetic field applied antiparallel to the anisotropy axis {6h = 7r) as 

2 13/2 x 3I/2 y 



r = 



7T 



1/2 



h 



K 1 S~ 1/2 e(l+4:K 2 ) 



1 - 16K 2 - 6AK 3 e - § 



(55) 



1 - 48^2 + 96 (K 3 - K'^j 

where the WKB exponent or the classical action S c i is shown in Eq. (52). Eq. (50) shows 
that in this case \<fi\ 1 is not valid, and therefore the problem can not be reduced to the 
one-dimensional motion problem. And the effective potential energy and the effective mass 
in one-dimensional form are not appropriate for the present case. 

Now we discuss the range of angles that Eq. (46) is valid. Introducing X = H — vr/2 
and 2 = 7T - Oh, from Eqs. (39), (46) and (52), we find that 1 ss (5 6 x 2~ 21 / 2 x 3~ 15 / 2 ) e 3 / 2 
and 02 ~ (5~ 6 x 2 39 / 2 x 3 15 / 2 ) e 3 / 2 . This means that Eq. (46) is almost valid in a wide range 
of angles 91° <0 H < 179° for e = 0.001. 
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For the single-domain FM nanoparticle with hexagonal crystal symmetry in the presence 
of an external magnetic field at arbitrarily directed angle, by using Eqs. (39) and (43) for 
H = tt/2, Eqs. (46) and (49) for tt/2 + O (e 3 / 2 ) < H < tt - O (e 3 / 2 ), and Eqs. (52) and 
(55) for H = tt, we obtain the ground-state tunnel splitting for MQC and the tunneling 
rate for MQT of the magnetization vector. Our results show that the tunnel splitting and 
the tunneling rate depend on the orientation of the external magnetic field distinctly. When 
Oh = tt/2, the magnetic field is applied perpendicular to the anisotropy axis, and when 
Oh = 7r, the field is antiparallel to the anisotropy axis. It is found that even a very small 
misalignment of the field with the above two orientations can completely change the results 
of tunneling rates. Another interesting observation concerns the dependence of the WKB 
exponent or the classical action with the strength of the external magnetic field. In a wide 
range of angles, the e (= 1 — H/H^j dependence of the WKB exponent S c i is given by e 5//4 , 
not e 3 / 2 for H = 7r/2, and e for H = 7r. Therefore, both the orientation and the strength 
of the external magnetic field are the control parameters for the experimental test for MQT 
and MQC of the magnetization vector in single-domain FM nanoparticles. 



V. CONCLUSIONS 



In summary we have investigated the tunneling behaviors of the magnetization vector 
in single-domain FM nanoparticles in the presence of an external magnetic field at arbi- 
trarily directed angle. We consider the magnetocrystalline anisotropy with the trigonal 
crystal symmetry and that with the hexagonal crystal symmetry. By applying the instanton 
technique in the spin-coherent-state path-integral representation, we obtain both the WKB 
exponent and the preexponential factors in the tunnel splitting between energetically degen- 
erate states in MQC and the tunneling rate escaping from a metastable state in MQT of the 
magnetization vector in the low barrier limit for the external magnetic field perpendicular to 
the easy axis (0 H = tt/2), for the field antiparallel to the initial easy axis {Oh = ?r), and for 
the field at an angle between these two orientations {tc/2 + O (e 3//2 ) < H < tt — O (e 3 / 2 )). 
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One important conclusion is that the tunneling rate and the tunnel splitting depend on the 
orientation of the external magnetic field distinctly. Another interesting conclusion concerns 
the field strength dependence of the WKB exponent or the classical action. We have found 
that in a wide range of angles, the e (= 1 — H/H c ) dependence of the WKB exponent or 
the classical action S c i is given by e 5//4 , not e 3//2 for Oh = vr/2, and e for 9h = 7r. We have 
obtained the temperatures corresponding to the crossover from quantum to thermal regime 
which are found to depend on the orientation of the external magnetic field distinctly. As a 
result, we conclude that both the orientation and the strength of the external magnetic field 
are the controllable parameters for the experimental test of the phenomena of macroscopic 
quantum tunneling and coherence of the magnetization vector in single-domain FM nanopar- 
ticles with trigonal and hexagonal symmetries at a temperature well bellow the crossover 
temperature. We have analyzed the validity of the semiclassical approximation performed 
in the present work, and have found that the semiclassical approximation should be already 
rather good for the typical values of parameters for single-domain FM nanoparticles. 

Recently, Wernsdorfer and co-workers performed the switching field measurements on 
individual ferrimagnetic and insulating BaFeCoTiO nanoparticles containing about 10 5 -10 6 
spins at very low temperatures (0.1-6K).! They found that above 0.4K, the magnetization 
reversal of these particles is unambiguously described by the Neel-Brown theory of thermal 
activated rotation of the particle's moment over a well defined anisotropy energy barrier. 
Below 0.4K, strong deviations from this model are evidenced which are quantitatively in 
agreement with the predictions of the MQT theory without dissipation. I The BaFeCoTiO 
nanoparticles have a strong uniaxial magnetocrystalline anisotropyl However, the theoreti- 
cal results presented here may be useful for checking the general theory in a wide range of 
systems, with more general symmetries. The experimental procedures on single-domain FM 
nanoparticles of Barium ferrite with uniaxial symmetryB may be applied to the systems with 
more general symmetries. Note that the inverse of the WKB exponent B~ l is the magnetic 
viscosity S at the quantum-tunneling-dominated regime T <^ T c studied by magnetic re- 
laxation measurement s.0 Therefore, the quantum tunneling of the magnetization should be 
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checked at any Oh by magnetic relaxation measurements. Over the past years a lot of exper- 
imental and theoretical works were performed on the spin tunneling in molecular Mn 12 -Ac!! 
and FeJ^ clusters having a collective spin state S — 10 (in this paper S = 10 6 ). Further 
experiments should focus on the level quantization of collective spin states of 5* = 10 2 -10 4 . 
We hope that the theoretical results presented in this paper may stimulate more experiments 
whose aim is observing macroscopic quantum phenomena in nanometer-scale single-domain 
ferromagnets. 
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APPENDIX A: EVALUATION OF THE PREEXPONENTIAL FACTORS IN 

WKB TUNNELING RATE 

In this appendix, we review briefly the procedure on how to calculate the preexponential 
factors in the WKB rate of quantum tunneling of the magnetization vector in single-domain 
FM particles, based on the instanton technique in the spin-coherent-state path-integral 
representation. @ The preexponential factors in tunneling rate (MQT) or the tunnel splitting 
(MQC) are due to the quantum fluctuations about the classical path, which can be evaluated 
by expanding the Euclidean action to second order in small fluctuations. Then we apply 
this approach to obtain the instanton's contribution to the ground-state tunnel splitting 
for resonant coherently quantum tunneling of the magnetization vector in FM particles 
with trigonal crystal symmetry in an external magnetic field applied perpendicular to the 
anisotropy axis (considered in Sec. Ill) in detail. 

In Ref. 12, Garg and Kim have studied the general formulas for evaluating both the 
WKB exponent and the preexponential factors in the tunneling rate or the tunnel splitting in 
the single-domain FM particles based on the instanton technique in the spin-coherent-state 
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path-integral representation, without assuming a specific form of the magnetocrystalline 
anisotropy and the external magnetic field. Here we explain briefly the basic idea of this 
calculation. Such a calculation consists of two major steps. The first step is to find the clas- 
sical, or least-action path (instanton) from the classical equations of motion, which gives the 
exponent or the classical action in the WKB tunneling rate. Instantons in one-dimensional 
field theory can be viewed as pseudoparticles with trajectories existing in the energy barrier, 
and are therefore responsible for quantum tunneling. The second step is to expand the Eu- 
clidean action to second order in the small fluctuations about the classical path, and then 
evaluate the Van Vleck determinant of resulting quadratic form. For single-domain FM 
particles, writing 9 (r) = 9 (r) + 6>i (t) and <p ( T ) — 4> ( r ) + </>i ( r ) ; where 9 and (f> denote the 
classical path, one obtains the Euclidean action of Eq. (2) as Se [9 (t) , (r)] ~ S c i + S 2 S 
with Sd being the classical action or the WKB exponent and 5 2 S being a functional of small 
fluctuations 9i and 



where E„ = (d 2 E / d9 2 ) e=M , E H = (d 2 E/d9dfa e=u= ^, and E H = (d 2 E/d<J> 2 ) e= - e ^. 



Under the condition that E^ > 0, the Gaussian integration can be performed over fa, and 
the remaining 9\ path integral can be casted into the standard form for a one-dimensional 
motion problem. As usual there exists a zero-mode, d9/dr, corresponding to a translation of 
the center of the instanton, and a negative eigenvalue in the MQT problem.0@ This leads to 
the imaginary part of the energy, which corresponds to the quantum escaping rate from the 
metastable state through the classically impenetrable barrier to a stable one. The resonant 
tunnel splittings of the ground state for the MQC problem can be evaluated by applying 
the similar technique. What is need for the calculation of the tunneling rate (in MQT) and 
the tunnel splitting (in MQC) is the asymptotic relation of the zero mode, d9/dr, for large 




(Al) 



T 



d9/dr w ae 



as ( 



oo. 



(A2) 
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The new time variable ( in Eq. (A2) is related to r as 

d( = dr/2A(e(T)^(r)), (A3) 

where 

A (0, 0) = HS 2 sin 2 9/2VE M . (A4) 

The partial derivatives are evaluated at the classical path. Then the instanton's contribution 
to the tunneling rate for MQT or the tunnel splitting for MQC of the magnetization vector in 
single-domain FM nanoparticles (without the contribution of the topological Wess-Zumino, 
or Berry phase term in the Euclidean action) is given 

|a|Gu/7r) 1/2 e-^. (A5) 

Therefore, all that is necessary is to differentiate the classical path (instanton) to obtain 
d9/dr, then convert from r to the new time variable ( according to Eqs. (A3) and (A4), 
and read off a and /i by comparison with Eq. (A2). If the condition > is not satisfied, 
one can always perform the Gaussian integration over 9\ and end up with a one-dimensional 
path integral over (f>\. 

Now we apply this approach to the problem of resonant coherently quantum tunneling of 
the magnetization vector between energetically degenerate easy directions in single-domain 
FM nanoparticle with trigonal crystal symmetry in an external magnetic field applied per- 
pendicular to the anisotropy axis. After some algebra, we find that 

E H « 2K X (1 + 12K 2 - e) , (A6) 

which is positive. So we can perform the Gaussian integration over <pi directly. The relation 
between r and the new imaginary-time variable ( for this MQC problem is found to be 

r = , ^ rC- (A7) 

2KxV (1 + 12K 2 - e) 

It is easy to differentiate the instanton solution to obtain 
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d5 




(A8) 



as C 



oo. Thus, 



a = 



(A9) 



and 



// = v^S f 1 - ^ 2 + J 



(A10) 



Substituting Eqs. (A9) and (A10) into the general formula (A5), and using Eq. (17) for the 
classical action or the WKB exponent, we obtain the instanton's contribution to the tunnel 
splitting HA as expressed in Eq. (21) for nanometer-scale single-domain ferromagnets with 
trigonal crystal symmetry in the presence of an external magnetic field applied perpendicular 
to the anisotropy axis. 

The calculations of the tunnel splitting and the tunneling rate of the magnetization 
vector for other MQT and MQC problems considered in the present work can be performed 
by applying the similar techniques, and we will not discuss them in any further. 
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Figure Captions: 

Fig. 1 The 8 (= 6 - 0„) dependence of the effective potential ~E X (5) for 9 H = n/2 (MQC). 

Fig. 2 The <5 (= — O ) dependence of the effective potential Ei(S) for = 37r/4 
(MQT). Here, K 2 = 0.001. 

Fig. 3 The 9h dependence of the relative classical action S c i (9 H ) / S c i (9h = 3n/A) in the 
trigonal symmetry with e = O.OOf and K 2 = O.OOf by numerical and analytical calculations. 
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